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1 INTRODUCTION 



At a recent post-seminar gathering, Herb Wilf casually mentioned to those of 
us assembled the fact that the quantities ((k) — 1 sum to 1, i.e., 



He then declared that when a sequence of nonzero positive numbers sums to 
1, the entries of the sequence should well be interpretable as the probabilities 
of something, and asked what it might be in this case. The challenge here is 
to provide events as interesting as the numbers £(fc) — 1 themselves: rich with 
internal relations and easy to describe. One criterion for a really good answer 
to this challenge would be that the response supplies us with an intuitive feel 
for the relations among the quantities ((h) — 1. 

The context of the discussion was number-theoretic, so after some consider- 
ation I formalized the challenge thus: 

Problem 1. Partition N into sets {^2,^3, • ■ ■ } such that the asymptotic den- 
sities lim n ^oo tt# (A k n {1, 2, . . . , n}) = ((h) - 1. 

This note provides one such partition. There are, it will be seen, others, but 
aesthetics guides the choice presented here: the partition can be described in 
terms of components that build on each other, and we can produce an algorithm 
that can quickly determine the Ak to which any number belongs. 




2 THE CLASSES B k 



The key both to constructing the classes Ak and to proving that they have the 
required densities is to break down the £(fc) — 1 thus: 



where the entries of the last line are the summations of the (infinite) columns 
above. The strategy is now to construct classes each of which have asymptotic 
densities given by the individual entries of the table. As those entries are all 
unit fractions, our task is much simplified: ask for a set of integers of some 
peculiar irrational density and intuition would not necessarily leap to service; 
but ask for one-sixth of the integers, and any residue class mod 6 will do. So 
we begin in this section with that last row, by producing classes Bj of densities 
5' 5> 11 > • • • ' re P resen thig each column. In the next section we break each Bj 
apart into sections Bj,k of densities p-, 2 < j, k < oo, to produce the individual 
entries of the table. The finishing stroke is then to collect across rows to form 
the Ak- 

The constructions we use are often explained nicely by recasting them in 
terms of the factorial-base representation of a number, its expansion as x — 
di • l! + <22 • 2! + • — h a n n\, < a t < i). This expansion is unique, as the student 
unfamiliar with this base system is encouraged to try proving on their own. For 
more on such expansions, see pQ. 

The class B2 must have asymptotic density \. Let us select it to consist of 
those x = 1 mod 2. 

Modulo 6, the classes 1,3, and 5 have been assigned, so for B3 we select the 
residue class 2 mod 6. 

Modulo 24, we require two residue classes in order that B4 should have 
density j^. We select the first two available: those x = 4 or 6 mod 24. 

Modulo 120, we require six residue classes in order that B5 should have 
density and the first six available are 10, 12, 16, 18, 22, and 24. 

Notice that we cannot use modulus 20 for S5: every residue class mod 20 has 
already had some subprogressions assigned. Instead, we will use the moduli jl. 
Since for m < j, residue classes modulo ml are the same as ^ equally spaced 
residue classes modulo jl, all previous assignments can be viewed as multiple 
assignments modulo jl. For each Bj, we select the first (j — 2)1 residue classes 
that have not yet been assigned. What are these in general? 

Lemma 1. The assigned residue classes for B m are 

m-2 

{x = 0! + aj ■ jl raodral , 1 < a j < j} ■ 

3=1 

Proof. In the examples above, by the time we have selected component residue 
classes for Bj we have assigned membership in some Bi to the natural numbers 



1, 2, . . . , (j — 1)!, as part of residue classes either modulo jl or modulo divisors 
thereof (previous factorials). Let us proceed by induction, and suppose that this 
has been the case through B m -\: for each b < m—1, the numbers 1,2,... (6— 1)! 
have been assigned to some Bj, j < b, as part of residue classes modulo jl. 

Given this induction hypothesis, what is the first residue class modulo ml 
available for assignment to B m l None of 1 through [m — 2)!, to begin with: 
each of those classes were assigned to B m ~\ or some previous class. We do 
know that none of the next (to — 2)! residue classes mod ml, (to — 2)! + 1 
through (m — 2)! + (to — 2)!, are assigned to B m -\. Examine this segment. 

Since (to — 2)! = mod (to — 2)!, our assumption tells us that the numbers 
(m — 2)! + 1 through (to — 2)! + (to — 3)! must have been assigned to B m -2 or 
earlier classes, as residues modulo (to — 2)! or divisors thereof. However, none 
of (to — 2)! + (to — 3)! + 1 through (to — 2)! + (to — 3)! + (to — 3)! were assigned to 
£>m-2, and this segment is wholly contained within those numbers we already 
knew were not assigned to B m -i- At each stage of the analysis, we have an 
unbroken segment of assigned classes for the numbers 1, 2, . . . , (to — 2)! + (to — 
3) ! + • • • + (to — i) I , and we know that the next (to — i)l numbers are not assigned 
to Bi through B m —\. At the last step, we find that the last number assigned is 
(to — 2)! + (to — 3)! + • • • + 2! + l!, an odd number in class B 2 - 

The very next number is in the segment we know has not been assigned 
to #2 or any previous class, so it is available. Thus the first available residue 
class modulo ml to assign to B m is (to — 2)! + (to — 3)! + • • • + 2! + 1! + 0!, 
a construction for which there are unfortunately multiple labels and notations 
in current usage. The relevant sequence 1,2,4,10,34,154,874,... is Sloane's 
A003422 ([3]) missing a leading zero term, and is there called the left factorial 
and denoted In = 0! + • — h (n — 1)!. We will use this notation, so that the first 
unassigned class modulo ml is !(to — 1). 

Now let us determine the remainder of the assigned classes mod ml. At 
the next-to-last step of our analysis previously, assuming to > 4 so that such a 
step exists, we determined that ! (to — 1) — 2 must have been an assigned class 
modulo 3!, that is, ! (to — 1) — 2 = 2 mod 6. Then ! (to — 1) is not such a class 
and neither is ! (to — 1) + 2: these are the two even classes mod 6 that were 
not assigned to B3. Nor can they be assigned to Bi for i > 3 as part of residue 
classes modulo multiples of 3!. Thus, we can select these two as our first two 
residue classes modulo ml for B m . If to > 5, we can ascend backward another 
step: there we find that !(to — 1) — 4 and !(m — 1) — 6 had to be the assigned 
residue classes mod 4!. We already knew l(m — 1) and !(to — 1) + 2 were not in 
those classes, and here we find that they can be joined by ! (to — 1) + 6 or +12, 
and !(to — 1) + 2 + 6 or +12. These are the six residue classes modulo 24 that 
were left after we assigned two of them for B4. 

At the i-th step of the ascent, we find that to all the entries we previously 
picked, we can add as many as i — 1 multiples of il. Thinking about it in the 
factorial-base representation, using 010203 ... to denote ai-l!+02-2!+03-3! + . . . , 
the assigned classes will be 02111 • ■ • + O62O364&5 . . . , with < bi < i — 1. And 
these are exactly the numbers claimed in the Lemma. □ 
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2.1 A Fast Membership Algorithm 

Having assigned the whole numbers to their various B m , we would like a means 
of determining which B m a given x belongs to without constructing every B m 
until we assign x. As it turns out, systematically selecting the first available 
classes at each step helps us write a computationally simple algorithm to do so. 

Write x in the factorial-base representation, x — x% ■ 1! + £2 • 2! + X3 ■ 3! + 
• • • + Xj ■ jl, < x% < i. We know that if the least positive residue of x modulo 
any n\ is less than (n — 1)!, x must have been assigned to one of the Bj< n , 
because from our proof above, the assigned residue classes for the Bj< n were an 
unbroken string from 1 to (n — 1)!. 

A " carry" in factorial-base addition occurs when the sum of the coefficients 
on the two summands on il is at least i + 1, since (i + l)il — 1 ■ (i + 1)1. But 
the Lemma tells us that the residues assigned to B m are of the form 02111 ••• + 
• ■ ■ , with < bi < i — 1. So in order for a to ever occur at the il place 
of the factorial-base expansion of a residue assigned to B m , when 1 < i < m— 1, 
we need a very specific summand. We need 62 = 1, so that we carry from 2!, 
leaving a 0. We then need 63 = 2, so that the sum of coefficients on 3! is 1 from 
!(rn — 1), plus 2 from 63, plus a carried 1, making 4, so we carry 1 and leave a 
again. This string has to continue up to the il term where we desire a 0. 

Thus we have one of two cases that will diagnose where !(m — 1) ends for a 
given x. We have that the least positive residue of x modulo ml is < (m — 1)!. 
Possibly x = (m — 1)! modulo ml exactly, so the string of x.- L starts with a 
possibly empty string of 0s that terminates with a 1 in position x m -\\ we added 
the largest possible value for every bi, and carried at every step of the addition. 
If we did not carry, then there is some smallest m such that x m _i = but no 
previous Xi = 0, except for possibly an initial string that does not terminate 
with a 1. Write out x in the factorial-base representation, check for the first 
behavior, and if it doesn't happen find the first 0. The index where whichever 
of these occurs, occurs, is m — 1. 

3 THE CLASSES A k 

Now that we have identified the classes B m , there remain the tasks of breaking 
them up into subclasses to give us the individual entries of TableQ] and collecting 
the corresponding parts "horizontally" to form the At- 

Each of the columns' individual entries decrease in geometric progression. 
The class Bi has density \ — ^ + j + ^ + ■ ■ ■) ■ The class B3 has density 

\ = \ (I + I + J? + • ••)• The class Bi has densit y T2 = 15 (I + Ts + • • 0' 
etc. 

These geometric series give us the subclasses we need. We assign the earliest 
fraction m=± of B m to the first subset, the next fraction — ^ m ~ 1 - ) to the next set, 

171 m ' mm 1 

and so forth. This can be done by a one-step digit test. For an example from 
familiar territory, to obtain ^ of the integers, take those ending in 1 through 9. 
Then to obtain of the integers by choosing yjr of the remaining tenth, take 
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those that end in but have one of 1 through 9 in the tens place. 

Since B m is made up of residue classes modulo ml, rewrite x as x = xq + 
ml ■ (b^mP + b\m x + b2m 2 + . . . J , < hi < m — 1. If n is the smallest number 
such that b n ^ m — 1 assign x to the set B m . n +2- hi table [U x is part of the 
set with asymptotic density given by the entry in column m, row n + 1, which 
represents an entry in the sum for A n +2- 

Thus, for example, we break up B2 as follows: 82,2 = {x = 1 mod4(= 2-2!)}, 
with density j, 82,3 — {x = 3mod8(= 4 • 2!)}, with density i; 6 2 ,4 = {x = 
7mod 16(= 4 • 2!)}, etc. We break up B3 into fractions of size |, |, etc.: B% : 2 — 
{x = 2or8modl8(= 3 • 3!)}, B 3 ,3 = {x = 14or32mod54(= 9-3!)}, etc. 

We now sum up by collecting these subsets across rows: Ak = B m ,k- Since 
the B m ^k are disjoint and the series of partial sums of their densities converges 
absolutely to £(fc) — 1, Ak has exactly the required asymptotic density. 

The algorithm to determine which class Ak a given x belongs to is thus: 

Algorithm. Write x = ax ■ l! + 02 ■ 2! + 0,3 ■ 3! + . . . , < a» < i. 

Let j be the smallest index such that aj ^ 0. If aj = 1, then m = j + 1. If 
Cij =/= 1, then m is the smallest number such that aj 7^ for j < i < m — 1, and 
a m _i = 0. Then x £ B m . 

Rewrite x as x = Xq + ml ■ (bom + b\m x + fe 2 m2 + •••)> < b\ < m, 
< xo < ml. Let n be the smallest index such that b n ^ m — 1. Then k = n + 2, 
and x 6 Ak ■ 



4 DATA AND SPECULATION 



The first few of the sets Bk and Ak are listed below. 




B 2 = 


{1,3,5,7,9,11,13,15,...}. 




£3 = 


{2,8,14,20,26,32,...}. 




Bi = 


{4,6, 28,30, 52,54, 76,78,...}. 




B 5 = 


{10, 12, 16, 18, 22, 24, 130, 132, 136, 138, 142, 144, . . . }. 




B 6 = 


{34, 36, 40, 42, 46, 48, 58, 60, 64, 66, 70, 72, ... , 120, . . . }. 




A 2 = 


{1, 2, 4, 5, 6, 8, 9, 10, 12, 13, 16, 17, 18, 20, 21, 22, 24, 25, 26, . . 


■}■ 


Az = 


{3, 11, 14, 19, 27, 32, 35, 43, 51, 59, 67, 68, 75, 76, 78, 83, 86, . . 


■}■ 


M = 


{7, 23, 39, 50, 55, 71, 87, 103, 104, 119, 135, 151, 167, 183, 199, 


...}. 


A 5 = 


{15, 47, 79, 111, 143, 158, 175, 207, 239, 271, 303, 320, 335, . . . 


}• 


More 


terms, for the series up to A10, have been submitted 


to Sloane's 



database [5] and should be available by the time of publication of this note. 
In terms of their component residue classes, the sets are: 



B 2 


= {x 


= lmod2}. 


B 3 


= {x 


= 2 mod 6}. 


Bi 


= {x 


= 4, 6 mod 24}. 


B 5 


= {x 


= 10,12,16,18,22,24modl20}. 


B 6 


= {x 


= 34, 36, 40, 42, 46, 48, 58, 60, 64, 66, 70, 72, ... , 120 mod 720}. 


A 2 


= {x 


= lmod4; 2, 8 mod 18; 4, 6, 28, 30, 52, 54 mod96 . . . }. 


A 3 


= {x 


= 3 mod 8; 14, 32 mod 54; 76, 78, 172, 174, 268, 270 mod384 . . . } 
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M = {x = 7 mod 16; 50, 104 mod 162; 364, 366, 1134 mod 1536 ... }. 

Some variations of this construction could be explored. Building Bj, we had 
to make a series of choices. We chose the odd numbers to form B2, but the even 
numbers would have worked just as well. We chose the residue class x = 2 mod 
6 for the class #3, but could just as easily have chosen x = 4 or 6. Choosing 
the first available residue class for each factorial modulus gave us classes B m 
that could be determined with a short algorithm, providing a tidy answer to the 
original problem. However, other choices lead to answers with different features. 

For example, suppose that at each step we choose the last available residue 
class modulo j! to construct Bj. Choose those x = 2 mod 2 for £> 2 - Of the three 
available odd residue classes mod 6, choose x = 5 mod 6 for B3. Use the classes 
x = 19, 21 mod 24 for B4. With such a decision procedure, some numbers are 
never assigned to any Bj at all! 

The "missed set" is 1,3, 7, 9, 13, 15, 25, 27, 31,33, •••= 1 + £ a;£!, < a* < i. 
All of the Bj will still have the same densities, summing to 1, so the missed set 
is "small" in that it is of asymptotic density 0. In fact, given any arithmetic 
progression X mod Y, some subprogression will be assigned to a Bj. (Show it!) 
On the other hand, at any given point the missed set may be rather large for 
some purposes: 1 jn of the numbers smaller than nl are permanently unassigned. 

These two assignment procedures are in some sense on opposite poles of an 
entire ensemble of possible procedures. An interested reader might burrow a 
layer deeper than we have: assign some straightforward process for describing 
and choosing assignment procedures to construct each Bj from residue classes 
mod j\ and examine the resulting ensemble of all possible constructions. Docs 
it possess any striking structural features? 

An early stab at this problem involved the powerfree numbers. Squarefree 
numbers are those whose prime factors are all distinct: they have density 
in the whole numbers. Cubefree numbers have factorizations in which no prime 
factor is repeated more than twice, so squarefree numbers are also cubefree. 
Cubefree numbers are of density ^y, and those which are cubefree but not 

squarefree are of density ^jy — ^y- Those that are 4th-powerfree but not 
cubefree are of density — ^y, and so on. Every whole number is (fc + l)st- 
powerfree but not fcth-powerfree for some fc, the largest exponent in its prime 
factorization. 

These densities seemed to suggest the possibility of defining a simple prob- 
ability distribution on N that assigned a total probability of £(fc) — 1 to the 
event that a random integer variable would be (fc + l)st-powerfree but not fcth- 
powerfree. It is easy to provide such a distribution by brute force - say, giving 
the nth number which is (fc + l)st-powerfree but not fcth-powerfree a probability 
of (C(fc) — 1) • ^rO But this does not seem to be a particularly illuminating il- 
lustration of any relations between the values £(fc); indeed, such a construction 
works with any partitioning of the integers, and any set of probabilities instead 
of £(fc) — 1. A distribution based on the factorization of x would seem much 

x Typo in this line corrected in preprint at kind communication from Michael Lugo. 
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more natural; can a simple one be produced? 

In closing, I would like to mention a personal recollection. In general, the 
density of Bj is j7jz?n ■ Seeing it written on the wall of my office, I was reminded 
of the first place I ever encountered Leibniz' summation of the triangular series: 
William Dunham's delightful popular-mathematics text, Journey Through Ge- 
nius [2]. I read this book in high school, and it motivated in considerable part 
my decision to pursue mathematics in college. The present note gives me an 
opportunity to thank Mr. Dunham sincerely for the service. 
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